home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 22 / CU Amiga Magazine's Super CD-ROM 22 (1998)(EMAP Images)(GB)[!][issue 1998-05].iso / PowerPC / Programming / PPCsiod / sources / math3.c < prev    next >
Encoding:
C/C++ Source or Header  |  1993-10-01  |  5.5 KB  |  199 lines
  1. /* Scheme In One Define.
  2.  
  3. The garbage collector, the name and other parts of this program are
  4.  
  5.  *                     COPYRIGHT (c) 1989 BY                              *
  6.  *      PARADIGM ASSOCIATES INCORPORATED, CAMBRIDGE, MASSACHUSETTS.       *
  7.  
  8. Conversion  to  full scheme standard, characters, vectors, ports, complex &
  9. rational numbers, and other major enhancments by
  10.  
  11.  *      Scaglione Ermanno, v. Pirinoli 16 IMPERIA P.M. 18100 ITALY        * 
  12.  
  13. Permission  to use, copy, modify, distribute and sell this software and its
  14. documentation  for  any purpose and without fee is hereby granted, provided
  15. that  the  above  copyright  notice appear in all copies and that both that
  16. copyright   notice   and   this  permission  notice  appear  in  supporting
  17. documentation,  and that the name of Paradigm Associates Inc not be used in
  18. advertising or publicity pertaining to distribution of the software without
  19. specific, written prior permission.
  20.  
  21. PARADIGM  DISCLAIMS  ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING
  22. ALL  IMPLIED  WARRANTIES  OF MERCHANTABILITY AND FITNESS, IN NO EVENT SHALL
  23. PARADIGM  BE  LIABLE  FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR
  24. ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
  25. IN  AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT
  26. OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
  27.  
  28. */
  29.  
  30. #include <stdio.h>
  31. #include <error.h>
  32. #include <string.h>
  33. #include <ctype.h>
  34. #include <time.h>
  35. #include <setjmp.h>
  36. #include <signal.h>
  37. #include <math.h>
  38. #include <limits.h>
  39.  
  40. #include "siod.h"
  41.  
  42. LISP plus(LISP args)
  43. {LISP l,s,res;
  44.  res=intcons(0);
  45.  for(l=args;NNULLP(l);l=cdr(l))
  46.    {s = car(l);
  47.     if (NNUMBERP(s))
  48.       err("+",s,ERR_GEN_ARG | ERR_NNUM);
  49.     res=plus2(res,s);}
  50.  return(res);}
  51.  
  52. LISP ltimes(LISP args)
  53. {LISP l,s,res;
  54.  res=intcons(1);
  55.  for(l=args;NNULLP(l);l=cdr(l))
  56.    {s = car(l);
  57.     if (NNUMBERP(s))
  58.       err("*",s,ERR_GEN_ARG | ERR_NNUM);
  59.     res=times2(res,s);}
  60.  return(res);}
  61.  
  62. LISP difference(LISP args)
  63. {LISP l,s,res;
  64.  if(NULLP(args))
  65.   err("-",NIL,ERR_FIRST | ERR_NNUM);
  66.  if(NULLP(cdr(args)))
  67.    res=flocons(0.);
  68.  else
  69.    {res=car(args);
  70.     if (NNUMBERP(res))
  71.       err("-",res,ERR_GEN_ARG | ERR_NNUM);
  72.     args=cdr(args);}
  73.  for(l=args;NNULLP(l);l=cdr(l))
  74.    {s = car(l);
  75.     if (NNUMBERP(s))
  76.       err("-",s,ERR_GEN_ARG | ERR_NNUM);
  77.     res=minus2(res,s);}
  78.  return(res);}
  79.  
  80. LISP quotient(LISP args)
  81. {LISP l,s,res;
  82.  if(NULLP(args))
  83.    err("/",NIL,ERR_FIRST | ERR_NNUM);
  84.  if(NULLP(cdr(args)))
  85.    res=flocons(1.);
  86.  else
  87.   {res=car(args);
  88.    args=cdr(args);}
  89.  for(l=args;NNULLP(l);l=cdr(l))
  90.    {s = car(l);
  91.     if (NNUMBERP(s))
  92.       err("/",s,ERR_GEN_ARG | ERR_NNUM);
  93.     res=divide2(res,s);}
  94.  return(res);}
  95.  
  96. LISP remainder(LISP x,LISP y)
  97. {LISP res;
  98.  double tmp;
  99.  x = tofloat(x);
  100.  y = tofloat(y);
  101.  if (NFLONUMP(x)||(modf(FLONM(x),&tmp)!=0.)) 
  102.   err("remainder",x,ERR_FIRST | ERR_NINT);
  103.  if (NFLONUMP(y)||(modf(FLONM(y),&tmp)!=0.)) 
  104.   err("remainder",y,ERR_SECOND | ERR_NINT);
  105.  res=flocons(fmod(FLONM(x),FLONM(y)));
  106.  return (res);}
  107.  
  108. LISP modulo(LISP x,LISP y)
  109. {LISP z;
  110.  double res;
  111.  x = tofloat(x);
  112.  y = tofloat(y);
  113.  if (NFLONUMP(x)||(modf(FLONM(x),&res)!=0.)) 
  114.   err("modulo",x,ERR_FIRST | ERR_NINT);
  115.  if (NFLONUMP(y)||(modf(FLONM(y),&res)!=0.)) 
  116.   err("modulo",y,ERR_SECOND | ERR_NINT);
  117.  res=fmod(fabs(FLONM(x)),fabs(FLONM(y)));
  118.  if(FLONM(y)<0)
  119.    res=-res;
  120.  z=flocons(res);
  121.  return(z);}
  122.  
  123. LISP expt(LISP x,LISP y)
  124. {LISP z;
  125.  errno=0;
  126.  if(RATNUMP(x))
  127.    {if(RATNUMP(y))
  128.       {if((RATNUM(x)<0))
  129.          {z = demoivre(pow((double)((double)RATNUM(x)/(double)RATDEN(y)),
  130.                            (double)RATNUM(y)),
  131.                        0.,
  132.                        (double)(1./(double)RATDEN(y)));}
  133.        else
  134.          z = ratcons(pow((double)RATNUM(x),
  135.                          (double)RATNUM(y)/(double)RATDEN(y)),
  136.                      pow((double)RATDEN(x),
  137.                          (double)RATNUM(y)/(double)RATDEN(y)));}
  138.     else
  139.       {y = tofloat(y);
  140.        if(FLONUMP(y))
  141.          {z = ratcons(pow((double)RATNUM(x),FLONM(y)),
  142.                       pow((double)RATDEN(x),FLONM(y)));}
  143.        else
  144.          err("expt",y,ERR_SECOND | ERR_NNUM);}}
  145.  else
  146.    {x = tofloat(x);
  147.     if(FLONUMP(x))
  148.       {if(RATNUMP(y))
  149.          {if((FLONM(x)<0))
  150.              {z = demoivre(pow(FLONM(x),(double)RATNUM(y)),
  151.                            0.,
  152.                            (double)(1./(double)RATDEN(y)));}
  153.           else
  154.               z = flocons(pow(FLONM(x),
  155.                           (double)RATNUM(y)/(double)RATDEN(y)));}
  156.        else
  157.          {y = tofloat(y);
  158.           if(FLONUMP(y))
  159.             {z = flocons(pow(FLONM(x),FLONM(y)));}
  160.           else
  161.             err("expt",y,ERR_SECOND | ERR_NNUM);}}
  162.     else
  163.       {y = tofloat(y);
  164.        if(NFLONUMP(y))err("expt",y,ERR_SECOND | ERR_NNUM);
  165.        if(COMPNUMP(x))
  166.          {z = demoivre((double)COMPRE(x),(double)COMPIM(x),FLONM(y));}
  167.        else
  168.           err("expt",y,ERR_FIRST | ERR_NNUM);}}
  169.  if(errno!=0)
  170.    raise(SIGFPE);
  171.  return(z);}
  172.  
  173. LISP demoivre(double ar,double ai,double e)
  174. {double tmp1,tmp2,tmp3,tmp4;
  175.  tmp1 = sqrt(pow((( ar*ar)+(ai*ai)),e));
  176.  tmp2 = atan2(ai,ar)*e;
  177.  tmp3 = cos(tmp2);
  178.  tmp4 = sin(tmp2);
  179.  return(compcons((float)(tmp1*tmp3),(float)(tmp1*tmp4)));}
  180.  
  181. LISP random(LISP x)
  182. {LISP z;
  183.  double r;
  184.  x = tofloat(x);
  185.  if (NFLONUMP(x)) err("random",x,ERR_GEN_ARG | ERR_NNUM);
  186.  r = (double) rand();
  187.  z = flocons(fmod(r,FLONM(x)));  
  188.  return (z);}
  189.  
  190. LISP randomize(LISP x)
  191. {
  192.  x = tofloat(x);
  193.  if (NFLONUMP(x)) err("randomize",x,ERR_GEN_ARG | ERR_NNUM);
  194.  if(FLONM(x)==0.)
  195.    srand(clock());
  196.  else
  197.    srand((unsigned long)FLONM(x));
  198.  return (truth);}
  199.